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Abstract 


General relativistic effects in the form of metric perturbations are usually 
neglected in the preheating era that follows inflation. We argue that in real¬ 
istic multi-field models these effects are in fact crucial, and the fully coupled 
system of metric and quantum field fluctuations needs to be considered. Met¬ 
ric perturbations are resonantly amplified, breaking the scale-invariance of the 
primordial spectrum, and in turn stimulate scalar field resonances via grav¬ 
itational rescattering. This non-gravitationally dominated nonlinear growth 
of gravitational fluctuations may have significant effects on the Doppler peaks 
in the cosmic background radiation, primordial black hole formation, gravi¬ 
tational waves and nonthermal symmetry restoration. 


98.80.Cq 


hep-ph/9808404 


Typeset using REVTgX 


1 






As recent ideas about the end of an inflationary era have shown, post-inflation reheating 
was one of the most violent and explosive processes occurring in the early universe H- 
The nonequilibrium, nonperturbative, resonant decay of the inflaton was demonstrated in 
Minkowski spacetime. The hrst steps toward a gravitationally self-consistent treatment, 
which modeled inflaton decay in an expanding, dynamical background spacetime, revealed 
that “preheating” may proceed with qualitative, as well as quantitative, differences from 
the Minkowski case |@J§]. Yet even these studies neglected an essential feature of gravita¬ 
tional physics: the production and amplihcation of metric perturbations attending such a 
sudden transfer of energy from the oscillating inflaton to higher-momentum particles. Other 
papers began the study of metric perturbations induced during preheating: in ^ a perfect- 
fluid analysis with Born decay was used, and crucial features of preheating were thus not 
incorporated; in these limitations are avoided, but the effects of the amplihed metric 
perturbations on the process of preheating itself are not considered. 

In this Letter, we pursue a more self-consistent relativistic treatment of preheating, by 
studying both the held huctuations (responsible for particle production) and the coupled 
metric perturbations (describing gravitational huctuations in the curvature). In the process, 
we clarify the question of causality and the amplihcation of long-wavelength perturbations 
during preheating. We give qualitative arguments, conhrmed by numerical results (see also 
0), to show that metric perturbations typically undergo rapid growth during multi-field 
preheating, and in turn act as a source and a pump for the growth of field fluctuations via 
gravitational rescattering, ultimately making the preheating process even more efficient than 
previously realized. Neglecting this general relativistic ehect can produce misleading results. 

Moreover, nonlinear growth of metric perturbations is itself of potentially major sig- 
nihcance, since it precipitates nonlinear density huctuations, leading to strong mode-mode 
coupling ehects, and nonlinear deviations from the conformally hat background. This could 
ahect observable quantities such as the cosmic microwave background (CMB) spectrum, 
complicating the usual predictions from inhationary models. Gravitational wave power 
could be enhanced by gravitational bremstrahlung, and primordial black hole production 
could occur without the need for special properties in the power spectrum. 

We work with the gauge-invariant formalism of to study the evolution of scalar per¬ 
turbations of the metric. (The amplihcation of gravitational waves at preheating has been 
considered in [Q.) For scalar perturbations, in the case of a scalar-held energy-momentum 
tensor and spatially hat background, the perturbed metric in the longitudinal gauge is 
= a?{ri) [(1 -|- 2<I)) drfi — (1 — 2$) dx^], where <F is the gauge-invariant gravitational po¬ 
tential. 

Realistic reheating models involve at least one held coupled to the inhaton (^(r;). First we 
consider models with only an inhaton, in order to give a simple illustration of some ehects. 
The gauge-invariant held huctuations 5(p and the metric perturbations <h obey the coupled 
Eqs. (6.40-42) of |^. Dehning the rescaled helds Y = a<I> and X = a5(p, these become, in 
momentum space: 



n'+ 



K + 


Xu = 2^"Yu , 




Yl = 

where V{‘p) is the potential. Id = a'/a and = SvrMpj^. 


( 1 ) 

( 2 ) 
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Eq. is a constraint, showing that if there is explosive growth in the held huctua- 
tions Xk - the heart of preheating - then this automatically constrains the gravitational 
huctuations to follow, and vice versa. This equation essentially says that if you “shake” 
the right-hand side of Einstein’s held equations then unavoidably you are 

simultaneously “shaking” the left hand side. Clearly, neglecting the metric perturbations 
Yk can be seriously misleading under many conditions, since it is tantamount to ignoring 
the perturbed Einstein equations. This may be reasonable in a slow-roll inhationary regime, 
but generally not in an oscillatory regime. 

A clear illustration is provided by super-Hubble modes in the simplest model, V = 
If 5(p^ is the held huctuation calculated by neglecting metric perturbations, then 


2 

for k 

x': + 


0, X* = a5(p* satishes Eq. 


with gravitational huctuations eliminated. 




X* = 0. This equation is then of precisely the same form as 


the background Klein-Gordon equation, so that 


oc (p ^ (posin( 6 r 7 ^)/( 6 p^), 


where and b are constants. The approximation arises from using the time-averaged 
scale factor d or. rf, and improves in accuracy as rj increases. When metric perturbations 
are incorporated, the long-wavelength solution is given in general by [see 0, Eq. (6.57)] 
6ip oc f a?dr]. We hnd that at reheating 


dip K. ^ipQ cos{br]^) , 

to lowest order in r]~^. Gravitational rescattering produces a non-decaying term in the held 
huctuations, which dominates the rapidly decaying huctuation 6ip^ calculated by neglecting 
gravitational huctuations. Relating in the usual way the held huctuations to the density 
of particles produced per mode, we hnd that this non-decaying solution indicates nonzero 
particle production. If gravitational perturbations are neglected, no particle production 
is found for this model. Gravitational rescattering dramatically alters the evolution of the 
matter-field fluctuations even in this simplest of all models, especially on super-Hubble-radius 
scales. Since preheating concerns primarily the behavior of such matter-held huctuations 
after inhation, it is thus crucial to study the coupled metric perturbation - held huctuation 
system. 

Now we consider the amplihcation of gravitational huctuations. The constrained system 
of equations (|l|)-( 0 ) has one degree of freedom, rehected in the decoupled equation Y^' — 
2 [ip" jpJ) Yl — l^K^ip''^ — Yk = 0. To avoid the periodic singularities when ip' = 0, Nambu 
and Taruya [^ employ the rescaled Mukhanov variable Q = a ( 2 -n)/(n+l) [X + [p>'/WjY], 
where n is the index in the power-law potential V = ipQ{ip/ipofi". Using the time- 
averaged forms of a and ip, they hnd that, to leading order in oc ^- 3 /(n+i)^ 


Q. 


1 -I- k'^Y 


T 


(4/r) sin2r Q = 0 , 


(4) 


where 7 is a constant. This equation has Mathieu form, y" + [Ak — 2qsm.2x\y = 0, with 
time-dependent Ak and q, so that modes can be drawn through instability bands by the 
expansion of the universe, if they were not already there initially i5|JS[]. Note that the 


1 -1 


resonance parameters scale as q/Ak or a ^2 _|_ 2)/(n+i) ^ where [3 is constant. 
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If metric fluctuations are neglected, then q/Ak oc a~^. Therefore, for n < 8, expansion is 
less effective in ending resonance when metric fluctuations are incorporated. 

Further discussion of single-held models is given in ||IO| , pd| , which conhrm the analytical 
conclusions of |^: in the n = 1 (quadratic potential) model, there is no resonance in the 
long-wavelength limit, whereas large, resonant growth is found for other models, such as the 
massless, quartically-coupled case (n = 2). However, as pointed out above, the single-held 
case is completely inadequate as a model of reheating, and we turn now to consider the 
multi-held case. 

The rapid growth of metric perturbations in the multi-held case (see Fig. 1) will produce 
a backreaction on the background quantities a and ip 


1^. Similarly, the amplihed held 
huctuations will grow to be of the same order as the tree-level terms, such as (p^, and hence 
will damp the inhaton’s oscillations, ending the parametric resonance. Let us consider an 
ehective single-held model. We may then estimate these two time-scales, expecting the initial 
preheating phase to end at rj^nd = minlr^m, pf}, where (<I>^(pm)) = 1 and (5(p^(pf)) = (p^. 
In a resonance band, <I>fc = Fk{ri)e^'^'^, where F). oscillates and pk is the Floquet exponent. 
From Eq. (^, and working in the saddle-point approximation, we may estimate 

($2) = {271)-^ J ~ «V'(V)/(4/iLa.) > 

where is the maximum Floquet exponent. For single-inhaton models of chaotic inhation 
with polynomial potentials, the slow-roll conditions are violated, and the inhaton begins to 
oscillate at Pq = aMp\, with a ~ 0.3. In models with a massive inhaton, we may further 
approximate |(p'| ~ amMpi. The held huctuations saturate their (linear-theory) upper limit 
at {Sp‘^) ~ Combining these yields (<F^(?7m)) ~ IQn'^a^nfl /The specihc spectrum 
of huctuations, governed by the values of pfc, depends on details of the potential. For modes 
subject to a parametric resonance, pk^a^ ~ cem, whereas modes subject to a negative- 

Hll 


coupling instability may have pkn 


0{m) 


Thus hrst-order analysis reveals that 


gravitational backreaction will become relevant at around the same time as the backreaction 
of the nonlinearly-coupled held huctuations. 

Another important issue is to demonstrate how super-Hubble modes may be amplihed 
causally. In inhation, modes which had been deep within the Hubble radius during inhation 
become amplihed and stretched to super-Hubble scales 00 . Reheating was believed not 
to be able to ahect these super-Hubble scales. In preheating, however, the coherence of the 
inhaton condensate immediately after inhation does allow for super-Hubble dynamics. At 
hrst, such behavior might appear to violate causality. This is not the case; consider the 
following: (1) the held equations are relativistic, hence causality is automatically built into 
the solutions; (2) causality concerns space-like related events, and does not translate into 
direct constraints on individual modes in Fourier space [^; (3) explicit calculation of the 
unequal-time two-point correlation function reveals that no mass or energy is being trans¬ 
ported superluminally by these super-Hubble resonances. See [|] for this calculation and 
further discussion. Similar conclusions regarding the possibility for the causal amplihcation 
of super-Hubble modes at preheating have been reached in IJOT. The main point we wish to 


emphasize is that causality restricts the shape of the spectrum of amplihed modes, but not 
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directly the wavelengths that can be amph£ed.|^ 

Preheating in single-held models is typically restricted to the narrow-resonance regime, 
and we may expect much larger effects in models with multiple scalar fields coupled to the 
oscillating infiaton, because the resonance parameter, q, may be much greater than unity 
The dynamics of such coupled-oscillator systems are in general chaotic and lead 
to enhanced particle production 


The multi-field generalization of equations 


can 


be given as 


3if$ + 


{k/aY + 


$ = 


1 

2 


K 




{5(pi)" -F SH^Sipi)' -f {k‘^/a?)5(pi = 
4<h(pi - 2Vi^ - T,Vij5ipj , 


(5) 

( 6 ) 
(7) 


dropping the mode label k, using proper time t, and writing Vi = dV/d^pi. Eq. m shows 
that resonant amplification of field fluctuations is accompanied by similar behavior of grav¬ 
itational fluctuations. In the single-field case, any resonant growth occurs with the same 
characteristics for metric and field fluctuations, and stability bands are thus also the same. 
In the multi-field case, this simple relation is broken, and the stability band structure is 
much more complicated. The non-infiaton fields (pi {i > 1) grow rapidly under resonance, 
while the infiaton ipi is strongly damped. Thus metric fluctuations $ grow more quickly 
than any of the field fluctuations 6ipi due to the quadratic products (pflipi. In models which 
include a substantial broad resonance regime for the coupled field fluctuations, resonance 
parameters are much larger than in the single-field case, typically in the range q 10^-10^ 
rather than g < 0(1) |^, so the associated Floquet indices are much larger 

These qualitative remarks are illustrated in Fig. 1 and amplified in the extensive simu¬ 
lations in [^, arising from integrating Eqs. i)-0 for the potential V = ^m‘^ip‘1 + ^g(p1ipl, 
which describes decay of the massive infiaton ipi into the boson field ip 2 - The graph shows 
the resonant amplification and nonlinear growth of gravitational fluctuations on both super- 
and sub-hubble scales. Also clear is the validity of the Floquet index as a characterizer of 
growth at strong resonance in an expanding universe (the slopes in the resonance bands are 
very nearly constant). The expansion has the effect of pulling modes through the resonance 
bands leading to phases of explosive growth and quiescence. Note that even with the ex¬ 
pansion of the universe included, the Floquet index is a useful concept at broad resonance - 
the growth of is almost exactly exponential in the resonance bands, as is evident in Fig. 
(1). A further crucial point, with wide-ranging implications, is that the non-gravitationally 
dominated evolution of perturbations during preheating breaks the scale-invariance of the 
primordial spectrum, as is evident from comparing the evolution of the k = 0, 20 modes in 

Fig. (1) §. 


^ The convenient approximation sometimes used in preheating, that the distribution of amplified 
modes falls as a spike, 6{k—kresona.nce), violates causality, since this requires that the field fluctuation 
contain correlations on all length scales, even for modes with k/a ^ H. 
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The numerical simulations in |P reveal that inclusion of the coupled metric perturbations 
enhances the strength of the resonances in Sipi, compared with when $ is neglected. This 
can be understood by the fact that the gravitational held has negative specihc heat. For 
this reason, <F does not act as a ‘parasite’ or ‘competitor’ with the d^pi for the energy of the 
oscillating inhaton, but rather can serve as a source or pump for held huctuations growth. 

Preheating in the narrow resonance and broad resonance regimes shows another impor¬ 
tant qualitative diherence: broad-resonance preheating in an expanding universe proceeds 
stochastically 0. The phases of the amplihed held modes are virtually uncorrelated between 
each moment when the oscillating inhaton passes through zero. Both the held modes and 
the metric perturbations will have stochastic driving terms. It has been shown that such 
stochastic terms in general remove all stability bands, so that all modes k > 0 grow sub¬ 
ject to a parametric resonance, with, in general, larger characteristic exponents than in the 


simple periodic case iT^,pT 


Note that once the huctuations have become strongly nonlinear, the linearized pertur¬ 
bation equations are no longer valid, and mode-mode coupling, arising due to convolutions 
which are ignored in the linearized Eqs. ®- 0 . must be included 0. 

Nonlinear gravitational huctuations could produce various observable signatures, in par¬ 
ticular on the CMB. Firstly, if the preheating phase is followed by a second round of inhation¬ 
ary expansion, then sub-hubble amplihed modes would get stretched into observationally- 
interesting scales. Such double-inhation is typical for most realistic scenarios based on 

0 . 


supergravity or supersymmetry 


Secondly, acoustic Doppler peaks at i > 100 could 


also be ahected. These are a key prediction of most inhationary models. In contrast, defect 
models often predict no secondary peaks and a shift in the position of the hrst peak, because 
the metric perturbations are produced randomly, receiving out of phase “kicks” due to the 
mode-mode coupling inherent in the nonlinearity of defect models El- 

In the multi-held, broad-resonance case, stochastic amplihcation leads to nonlinear met¬ 
ric huctuations, producing mode-mode coupling which drives the unequal-time correlation 
function A — > 0 for \ri — ri\ > r]c, the coherence time of the system. As pc ^ 0 we get 6- 
correlated (white) noise [|^], memory of initial conditions (and hence coherence) is lost, and 
the held evolution mimics the active, incoherent evolution of defects. The crucial question 
is whether the nonlinear mode-mode coupling survives local interactions and persists up to 
nucleosynthesis and photon decoupling. If so, there could be signihcant limits placed on 
inhationary reheating by observed element abundances and small-scale CMB anisotropies. 
By smearing out the Doppler peaks, surviving nonlinear mode-mode coupling would greatly 
reduce the ehectiveness of small-angle CMB observations for diherentiating inhationary sce¬ 
narios from defects models. The implied nonlinearity would lead to hybrid CMB anisotropies 
- passive and typically inhationary on large angular scales, active and defect-like on smaller 
scales. 

Another possible observational consequence arises from gravitational bremsstrahlung in 
scattering of the large scalar perturbations 0. Combined with our analysis of causality 
above, the extra gravitational wave power due to the previously neglected metric huctua¬ 
tions implies that the probability of detection in instruments such as LIGO may be higher 
than previously estimated, at least for preheating following chaotic inhation. Nonlinear 
gravitational huctuations could also create the large density contrast necessary to produce 
primordial black holes (PBHs) 1^, without the need for a blue spectrum or a large peak 
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in the power spectrum at some k. PBH limits may then be able to constrain the nature of 
preheating and the associated amplihcation of metric perturbations. 

Finally, amplihed metric perturbations affect non-thermal symmetry restoration at pre¬ 
heating by altering the inflaton’s effective potential via addition of the terms (see Eq. 
(10.68) of [0]): AV = . The hrst term was originally considered in 

while the second is due to the direct (gravitational) coupling between bip and <h. It is 
of the same perturbative order as the hrst term, and is missed if gravitational huctuations 
are neglected in preheating. 
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FIGURES 
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FIG. 1. Metric perturbation evolution in 2-field preheating, with q = gip\{tQ)/m? = 8 x 10^. 
The main graph shows the k = 0 mode, the inset the k = 20 mode. Initially, mto = 100 (well into 
the small amplitude phase; (fi{to)/mto = 3 x 10“^) and ^k{to) = 10“®, ^k{to) = 0. The k = 0 
mode becomes nonlinear at mt ~ 150 after less than 10 inflaton oscillations (i.e. well before the 
end of preheating) and continues growing without bound. Two strong resonance bands are evident, 
with different Floquet indices (given by the slopes of the dotted lines), for 140 < mt < 160 and 
250 < mt < 300. 


10 




